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ABSTRACT : 

This  is  the  third  year-end  report  on  the  EHD  contract. 

Work  with  the  steam  injector  continued  during  the  past  year. 
The  definition  of  the  slip  parameter  for  unsteady  flow  has  been 
used  in  conjunction  with  velocity  profiles  measured  with  a  hot  wire 
anemometer  to  evaluate  the  potential  role  of  EHD  (or  EGD)  in  anemometry. 

A  direct  measurement  of  particle  size  (using  holographic  interferometry) 

-k 
indicated  that  there  are  no  particles  of  size  10  m  or  greater.  A 

computer  program  describing  the  EHD  flow  has  been  written. 


FOREWORD 

This  is  a  continuation  of  the  EHD  Research  program  started  on 
September  of  I968  at  the  Naval  Postgraduate  School,  Monterey,  California, 
This  work  was  sponsored  by  the  Naval  Air  Systems  Command  under  the 
technical  cognizance  of  Dr.  H.  R.  Rosenwasser. 

One  Master  thesis  was  generated  and  it  is  included  in  the  Appendix. 
A  paper  entitled  "Unsteady-Flow  Electrogasdynamics"  is  presently  being 
submitted  for  publication.  A  Navy  patent  has  been  applied  for  the 
spectral  anemometry  principle.  A  progress  paper  was  given  at  the 
Energy  Conversion  Colloquium  held  during  9-10  November  at  the  Naval 
Postgraduate  School. 

The  help  of  LCDR  Bohley  and  Mr.  Patrick  Hickey  is  gratefully 
acknowledged.   The  laser  holographic  measurements  were  made  possible 
through  the  assistance  and  cooperation  of  Professor  Collins  of  the 
Aeronautics  Department,  NPS . 


II 


TABLE  OF  CONTENTS 

Page 

I.  Introduction  ----------------.---.__  i 

II.  Hot-wire  Anemometer  Readings  -------------  l 

III.  Optimum  Frequency  Response  --------------  3 

IV.  Laser  Holography  -------------------  5 

V.  Computer  Programs  ------------------  5 

VI.  Conclusions  and  Recommendations  -----------  7 

References  ------------------------- 

Appendix  —  Spectral  Anemometry  with  an  Electro- 

gasdynamic  Probe  ---------------- 

Distribution  List  --------------------- 


III 


LIST  OF  FIGURES  (Main  Text) 

Figure  Page 

1  Composite  of  Hot -Wire  Anemometer  Data  9 

2  EGD  Generator  Schematic  10 

3  Optimum  Frequency  Response  as  a  Function 

of  the  Radius  11 

k         Laser  Hologram  of  the  Vitual  Image  12 

5  Enlarged  View  of  Test  Region  12 

6  Two -Dimensional  Output  of  Potential  Distrihution  I3 

7  Two -Dimensional  Output  of  the  Electric  Field  ik 


LIST  OF  SYMBOIS  (Main  Text) 
E,    =  Breakdown  field  strength 
1e  1  =  Absolute  value  of  the  electric  field 
R    =  Particle  radius 
S((o)  =  Slip  parameter  of  frequency  u) 
V(ua)  -   R^B  velocity  at  frequency  ou 
p,    =  Mobility  of  particles 
T    =  Equilibration  time  for  particles 
u)    -   Frequency 


IV 


I.   imRQDUCTION 

During  the  past  year  progress  was  made  basically  along  three  fronts: 
First,  the  anemometry  application  was  the  subject  of  an  MS  thesis 
(included  here  as  the  Appendix).  Also,  a  paper  tying-in  the  work  on 
anemometry  and  the  EHD  generator  was  written  and  is  presently  being 
submitted  for  publication.   Next,  an  attempt  was  made  at  measuring 
particle  size  directly  (i.e.,  by  laser  holography).  Last,  a  computer 
program  which  solves  the  combined  Laplace -Poisson  equations  inside 
the  generator  was  written  and  checked-out. 

This  report  includes  the  final  hot-wire  anemometer  frequency  profile 
of  the  region  being  investigated.  As  pointed  out  in  last  year's 
report  ,  these  measurements  were  made  along  the  centerline  in  the  wake 
of  the  cylinder  at  a  location  5i™i  downstream  of  the  exit  plane  of  the 
injector  nozzle.  A  review  of  the  optimum  frequency  response  concepts 
follows  the  hot  wire  data.  Also  included  are  a  discussion  of  the 
laser -holograms  and  a  brief  summary  of  the  computer  program. 

II.   HOT-WIRE  AMIMOMETER  READINGS 

The  unit  that  gave  the  final  results  reported  here  is  a  Thermo 
Systems,  Inc.  model  1051-2  constant  temperature  hot  wire,  with  a 
frequency  response  of  50  KHz. 

The  hot  wire  readings  were  to  be  used  as  a  standard  of  comparison 
for  the  EGD  data  but  because  the  hot  wire  calibration  is  not  appropriate 
for  two  phase  flows,  air  was  substituted  for  the  steam  in  the  cylinder  . 
The  air  flow  was  matched  so  as  to  give  the  same  momentum  as  the  steam. 


Ttxe   absence  of  droplets  extended  the  life  of  the  hot  wires  which  was 

very  marginal  to  begin  with  due  to  the  intense  vorticity  of  the  wake 

of  the  cylinder.  Another  complication  arose  from  the  fact  that  the 

air  flow  through  the  steam  nozzle  was  typically  as  a  different  temperature 

from  the  main  air  flowing  around  the  cylinder.  Several  attempts  were 

made  at  equalizing  the  two  temperatures .  The  most  reliable  data 

seemed  to  indicate  that  the  nozzle  flow  does  not  have  a  measurable 

effect  on  the  main  flow.  This  holds  for  distances  of  at  least  5nim 

away  from  the  nozzle  exit  plane.  The  resulting  profile  is  shown  in 

Figure  1.   It  represents  the  main  air  and  the  nozzle  air  (8  psig) 

but  it  is  indistinguishable  from  the  hot  wire  profile  for  the  main 

air  itself.  There  are  two  dominant  frequencies  apparent  from  this 

figure,  namely,  1.2  KHz  and  10  KHz.  The  first  one  represents  twice 

the  Strouhal  frequency  for  the  cylinder  and  is  appropriately  so 

because  the  measurements  were  taken  along  the  centerline  of  the  wake 

where  the  vortices,  shed  off  the  cylinder  on  alternating  sides,  add 

up.  The  Reynolds  number  based  on  cylinder  diameter  for  this  flow 

(55  m/sec)  is  8x10  so  that  the  boundary  layer  on  the  cylinder  is 

laminar.  The  second  important  frequency  apparent  from  Figure  1 

seems  to  be  associated  with  the  nozzle  that  injects  steam  into  the 

wake,  see  Figure  2.  The  10  KHz  peak  appeared  whether  or  not  there 

was  air  issuing  from  the  nozzle.   The  exit  diameter  of  this  nozzle  is 

about  2mm  which  is  some  ten  times  smaller  than  the  cylinder  diameter 

and  perhaps  appropriately  the  frequency  is  nearly  10  times  higher. 

It  was  observed  further  that  10  KHz  is  apparently  a  high  enough 

frequency  so  that  the  dissipation  of  the  vortices  is  high;  at  a  station 


lOinin  away  from  the  nozzle  exit  plane,  the  10  KHz  peak  disappears 
whereas  the  1.2  KHz  peak  remains  nearly  as  intense.  So  the  fluid 
viscosity  can  act  on  these  "lesser  whorls"  but  the  flow  is  still  to 
he  considered  pseudo-turbulent  since  the  vortices  at  the  Strouhal 
frequency  are  basically  two  dimensional. 

There  are  no  other  significant  contributions  to  the  vorticity 
because  the  flow  turbulence  is  carefully  removed  upstream  and  all 
measurements  are  made  away  from  the  wall  boundary  layers .  There  is 
one  uncertainty  that  remains,  however,  and  that  involves  whether  or 
not  the  presence  of  large  particles  (i.e.,  two-phase  flow)  shifts 
the  intense  part  of  the  profile  towards  the  lower  frequencies .  The 
larger  particles  can  absorb  some  of  the  eddy  energy  more  or  less 
continuously  at  this  may  cause  a  significant  shift.  This  could  be 
described  with  a  different  "effective"  viscosity  for  the  two-phase 
flow.   The  only  way  to  resolve  this  dilema  is  to  make  measurements 
with  either  a  laser  or  an  ultrasonic  anemometer.   The  laser  device 
relies  on  the  presence  of  large  particles  and  measures  their  Doppler 
shift.  The  accuracy  of  these  devices  is  considered  below  that  of  the 
hot  wire  at  the  present. 

III.   OPTIMUM  FEEQUENCY  RESPONSE 

It  has  been  pointed  out  that  in  the  anemometer  application  of 
the  EGD  probe  the  charged  particle  size  determines  the  frequency 
response  of  the  device.  This  is  clear  insofar  as  the  inertia  of  the 
particle  governs  how  fast  it  can  equilibrate  with  the  eddies. 


In  last  year's  report,  it  was  shown  that  the  slip  parameter  (S) 
for  unsteady  flow  could  be  written  as 

S      =  r-T-y  +  T   tU  (1) 


max 


VW    V 


where   to    =  frequency 

[i,         =  particle  mobility 

K    =  breakdown  field  strength 

V(cu)  =  RIVB  intensity  of  the  free  stream  velocity  at  the 
frequency  uj 

T    =  equilibration  time  for  the  particles 

Since  the  mass  of  the  particle  only  appears  in  t  ,  the  question  arises 
as  to  what  may  be  a  desired  criterion  for  the  frequency  response  as  a 
function  of  size  which  si  consistent  with  a  small  slip  parameter.  This 
criterion  may  be  arbitrary,  i.e.,  fixing  uot^  to  be  some  small  fraction. 

Figure  3  shows  a  plot  of  cu  vs  radius  for  values  of  cot  equal  to 

-2      -3 
10   and  10  .  This  choice  fixes  u)  without  regard  to  the  value  of 

-2 
the  first  term  in  Equation  1  which  in  itself  may  be  higher  than  10  ; 

however,  a  high  slip  due  to  a  high  mobility  is  undesirable  for  either 

the  generator  or  the  anemometer.  Note  further  that  if  V((o)  is  very 

small  i.e.,  close  to  zero,  then  the  slip  may  be  very  high  but  since 

there  is  no  significant  contribution  by  the  flow  at  this  frequency  the 

slip  value  is  unimportant. 

Since  the  slip  parameter  plays  a  significant  role  in  the  generator 

efficiency,  it  is  necessary  to  determine  the  spectral  distribution  of 

the  flow  inside  the  generator  and  calculate  the  slip  at  each  of  the 

dominant  frequencies. 


IV.  LASER  HOLOGRAJHY 

In  last  year's  Final  Report  ,  results  from  charge  distribution 
measurements  seem  to  indicate  the  existence  of  charged  droplets  of 

o 

size  of  about  100  [x.  m.   This  is  because  about  3  x  10  electrons  per 
pulse  had  been  measured  at  a  steam  pressure  of  13  psig.   The  resolution 

that  we  could  obtain  with  the  laser  apparatus  was  itself  about 

2  3 
100  ^,m  and  the  limiting  factor  '  was  the  graininess  of  the  hologram 

reproduced  object.  As  can  be  seen  from  Figirres  k   and  5?  no  particles 

can  be  detected  and  we  conclude  that  the  particles  are  of  a  size  much 

smaller  than  100  \xm.       Note  that  the  nozzle  exit  diameter  is  O.lU 

inches  or  2.5  mm. 

Other  work  in  last  year's  Final  Report  seem  to  indicate  that  sizes 

were  in  the  order  of  10  ^.m  and  less.   The  direct  measurement  of  these 

particles  is  a  the  limit  of  the  technique  even  with  a  refined  set-up. 

V.  COMPUTER  PROGRAM 

k 
As  shown  by  Minardi  the  potential  distribution  within  an  EGD 

channel  is  described  by  a  "combined  Laplace-Poisson"  equation.  The 

two  equations  are  required  by  the  two-phase  space  charge  nature  of 

the  EGD  environment.  Poisson's  equation  models  the  space-charge 

region,  while  Laplace's  equation  models  the  space-charge  free  region. 

Electrode  geometry  and  space-charge  effects  on  the  potential  distribution 

are  of  interest  in  the  design  of  EGD  devices.  The  solution  of  the 

governing  equations  will  show  these  effects.  However,  a  closed-form 


solution  is  often  very  difficult  to  find,  and  any  attempts  at  a  solution 
will  be  complicated  "by  the  varying  boundary  conditions  imposed  by  the 
electrodes . 

In  order  to  have  an  analytical  input  to  the  design  process,  a 
numerical  solution  method  (successive  over -relaxation)  was  programmed 
for  the  IBM  360  computer  system.  Finite  difference  equations,  derived 
from  a  Taylor-series  expansion  about  the  point  in  question,  were  used 
to  build  the  numerical  model.  The  program  has  the  capability  for 
solving  axisymmetric  and  two-dimensional  channel  flows.  The  program 
output  is  designed  to  give  both  numerical  and  graphical  representations 
of  the  solution. 

Of  importance  in  the  computer  modeling  of  the  EGD  device  is  the 

5 
fluid  dynamics  involved.  By  means  of  a  similar  program,  Minardi 

has  shown  that  for  particle  sizes  of  interest,  the  fluid  dynamics  may 

be  uncoupled  from  the  electrodynamics.  Thus  it  is  possible  to  limit 

modeling  of  the  dielectric  gas  flow  to  a  constant  rate  of  jet  expansion 

and  insure  conservation  of  charge  within  the  jet. 

Input  variables  to  the  program  typically  include  those  describing 

the  electrode  geometry,  dielectric  characteristics,  charged  particle 

characteristics,  flow  control,  variables  matrix  dimensions  and  the 

output  desired.  The  program  simulates  a  steady  state  picture  of  the 

entire  EGD  situation  by  the  superposition  of  two  separate  matrices. 

One  matrix  will  carry  all  potential  information,  while  the  other  matrix 

carries  both  charge  distribution  and  some  flow  control  information. 

The  initial  part  of  the  program  is  devoted  to  setting  up  these 

matrices. 


Once  this  has  been  accomplished,  the  program  calls  for  an  iterative 
calculation  to  determine  the  potential  distribution,  and  a  calculation 
for  the  |e I  -field  distribution.  Remainder  of  the  program  is  devoted 
to  output  of  these  results. 

The  program  was  applied  to  the  geometry  shown  in  Figure  2.  The 
dielectric  gas  used  was  steam  issuing  at  397  m/sec  from  the  nozzle. 
This  corresponds  to  a  reservoir  conditions  of  8  psig  and  2i|0°F.   This 

is  a  slightly  super-heated  steam.  Particle  charging  was  considered  to 

-7  2 
be  corona  limited,  giving  a  particle  mobility  of  1  x  10   ra  /v-sec. 

-12 
A  dielectric  constant  of  8.9  x  10    coul/v-m  is  used.  Attractor 

potential  was  3  KV,  and  the  collector  potential  was  50  volts  below 

ground.  Figure  6  is  the  plotted  output  of  the  potential  distribution, 

showing  lines  of  constant  potential.  Figure  7  is  the  plotted  output 

of  |e  I  -field  distribution,  showing  lines  of  constant  value  of  E. 

The  potential  has  been  normalized  with  respect  to  attractor  or  corona 

potential,  and  the  E-field  values  with  respect  to  the  breakdown  field 

of  3  X  10  v/m.  The  geometry  is  represented  by  the  cross-hatched 

points . 

VI.   CONCLUSIONS  AM)   RECOMMENDATIONS 

It  is  evident  from  the  results  shown  in  the  Appendix  that  we  Lrve 
been  able  to  achieve  "acceptable"  responses  from  the  EGD  probe  up 
to  1000  Hz.  This  acceptability  is  based  on  the  correctness  of  the 
hot-wire  data  used  for  comparison.  From  Figure  3?  this  implies  that 
particle  sizes  are  centered  about  1  [im  or  higher  radii.  From  the 
holographic  measurements,  we  know  the  particles  are  not  much  above 


10  |j,m.   Our  efforts  at  changing  the  superheat  of  the  steam  in  order 
to  affect  particle  size  have  indicated  that  this  is  not  the  only 
relevant  parameter;  apparently  what  happens  to  the  charged  particles 
outisde  of  the  nozzle  is  equally  important.  There,  of  course,  the 
behavior  is  somewhat  random  and  a  statistical  sample  is  needed. 
Lowering  the  pressiore  from  13  psig  to  8  psig  is  a  step  in  the  right 
direction  "but  apparently  something  happens  that  gives  multiple  peaks 
in  the  distribution  of  particles.  These  results  point  out  two  things; 
First,  better  humidity  control  and  some  dilution  of  the  steam  inside 
the  cylinder  (through,  say,  a  bleed  of  nitrogen)  may  be  needed  to 
achieve  a  reduction  of  particle  size.  Second,  the  present  response 
is  sufficient  to  measure  some  interesting  turbulent  flows  which  need 
not  be  as  complicated  as  the  wake  of  the  cylinder  (for  example, 
grid-generated  turbulence). 

The  computer  study  will  be  directed  toward  the  design  of  guard- 
electrodes  inside  the  generator  channel  which  will  improve  the  output 
of  the  device.  The  program  is  further  useful  in  pointing  out  regions 
of  high  electric  field.   This  is  quite  significant  for  any  EHD  power 
generator . 

The  definition  of  the  coupling  of  the  charged  particles  with  an 
unsteady  flow  field  remains  as  an  important  goal  of  this  project. 
Most  EHD  flows  are  unsteady  to  some  extent  and  a  low  slip  parameter 
for  steady  flow  may  not  be  good  enough  if  significant  velocity 
components  exist  at  non-zero  frequencies. 
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FIGURE  1.      COMPOSITE  OF  HOT-WIRE  ANEMOMETER  DATA 
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FIGURE  h     LASER  HOLOGRAM  OF  VIRTUAL  IMAGE 


FIGURE  5  ENLARGED  VIEW  OF  TEST  SECTION 
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ABSTRACT 

This  study  involves  the  measurement  of  RMS  velocities 
covering  the  frequency  spectrum  25  to  20,000  hertz  using 
an  electrogasdynamic  (EGD)  probe.   Previous  work  proved 
the  feasibility  of  this  investigation.   The  flow  of  air 
behind  a  circular  cylinder  was  measured.   The  freestream 
velocity  was  185  ft/sec.   Charged  water  droplets  entering 
the  wake  through  a  jet  in  the  cylinder  at  sonic  velocity 
acted  as  the  media  for  conveying  information  to  the  probe. 
Hot-wire  anemometer  readings  made  in  the  same  spectral 
range  were  used  as  the  standard  of  comparison.   However, 
these  need  further  refinement. 

Small  droplet  formation  was  optimized  by  determina- 
tion and  control  of  the  relative  humidity  of  the  freestream, 
steam  temperature  and  steam  pressure.   Optimum  droplet 
charge  was  then  determined  and  controlled  with  corona 
voltage.   A  10  to  15  decibel  increase  in  high  frequency 
response  and  stable  performance  resulted. 
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I.   INTRODUCTION 

Previous  work  in  Electro^^asdynamics  has  indicated  a 
definite  correlation  between  EGD  probe  readings  of  col- 
lector current  and  the  mean  and  rms  velocities  measured 
by  a  hot-wire  anemometer  in  the  turbulent  wake  of  a  cir- 
cular cylinder  [6].   This  thesis  is  an  attempt  to  improve 
the  spectral  correspondence  of  EGD  probe  measurements  to 
hot-wire  results.   The  probe  is  a  stainless  steel  rod  with 
a  sharpened  tip  which  is  placed  in  a  moving  stream  of  air 
where  it  collects  charged  particles.   For  charged  particles 
of  proper  size  and  charge,  sufficient  coupling  may  result 
in  the  turbulent  flow  to  allow  measurement  of  the  spectrum 
of  turbulence  intensity. 

The  advantages  of  the  EGD  probe  as  a  diagnostic  tool 
over  the  hot-wire  anemometer  are  manifold.   The  probe  is 
very  rugged,  allowing  probing  of  intense  turbulent  fields. 
It  is  very  simple  and  inexpensive  and,  therefore,  can  be 
easily  manufactured.   The  probe  tip  area  is  small  and  sym- 
metrical, allowing  better  spatial  resolution  than  the  hot- 
wire.  Since  the  very  small  charged  particles  can  react  to 
high  frequency  oscillations,  the  frequency  response  of  the 
EGD  probe  may  be  higher  than  that  of  the  hot-wire.   In  ad- 
dition, signals  from  the  EGD  probe,  which  are  electrical 
in  nature,  need  little  or  no  further  amplification  [7]. 
Effects  corresponding  to  the  thermal  inertia  of  hot-wires 
are  not  present.   There  is  no  requirement  to  feed  the 
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output  through  a  linearizer  since  the  EGD  output  is  already- 
linear,  i.e.,  a  direct  measurement.   Also  eliminated  is  the 
effect  that  free  convection  currents  have  on  the  hot-wire 
zero  setting.   In  the  case  of  two-phase  fluids,  the  EGD 
probe  needs  no  separate  calibration  and  will  withstand 
the  impact  of  charged  droplets  or  particles  whereas  the 
fragile  hot-wire  will  readily  break  under  these  conditions. 
The  main  disadvantages  of  the  EGD  mechanism  are  the  require- 
ments for  a  source  of  charged  particles  of  optimum  size  and 
charge  and  for  a  high  voltage  supply  (along  with  the  neces- 
sary safety  precautions) .   The  optima  appear  to  be  rather 
critical . 

The  specific  objectives  of  this  thesis  are:   1)   im-' 
provement  of  the  high  frequency  response  of  the  EGD  probe 
[5];  2)  the  elimination  or  explanation  of  inconsistencies 
between  the  EGD  results  and  the  hot-wire  results;  and  3) 
the  determination  of  the  effect  that  various  parameters, 
i.e.,  humidity,  steam  pressure  and  temperature,  the  ratio 
of  collector  current  to  total  current,  etc,  have  on  the  re- 
sults.  A  lesser  objective  was  finding  a  means  of  facili- 
tating the  reduction  of  data  for  comparison  to  hot-wire 
data. 

This  work  complements  part  of  a  comprehensive,  long- 
term  project  being  conducted  by  Assistant  Professor  Oscar 
Biblarz,  Department  of  Aeronautics,  Naval  Postgraduate 
School.   Existing  EGD  equipment  installed  in  Building  230 
at  the  Naval  Postgraduate  School  was  utilized. 
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II.   THEORY 

In  1959  Bennett  [2]  pointed  out  that  by  carrying  elec- 
trical charge  in  the  volume  of  a  stream  of  air  instead  of 
on  the  surface  of  a  belt,  a  gaseous  Van  de  Graff  generator 
could  be  built.   Since  a  greater  volume  of  charge  would  be 
carried  at  higher  speeds,  much  higher  current  should  re- 
sult.  This  marked  the  beginning  of  the  EGD  generator. 

The  basic  EGD  mechanism  involves  particles  which  are 
introduced  into  an  air  stream  or  other  suitable  gas  at  the 
entrance  of  the  test  section  (Figure  1) .   For  the  source 
of  charges  an  aerosol  flow  in  conjunction  with  an  ionizer 
are  used.   Saturated  steam  is  injected  through  a  corona 
discharge  into  the  "turbulent"  wake  of  a  cylinder.   The 
corona  discharge  is  produced  by  applying  a  high  voltage 
between  a  sharp  needle  and  an  attractor  ring  which  is  em- 
bedded in  the  wall  of  the  steam  nozzle  located  on  the 
downstream  side  of  a  right  circular  cylinder.   Electrons 
from  the  negative  needle  drift  toward  the  attractor  and 
attach  themselves  to  the  then-growing  water  droplets  form- 
ing negative  charged  particles.   Some  of  these  charged 
particles  are  in  turn  dragged  or  convected  by  the  air  flow 
down  the  conversion  channel  where  their  electrons  are 
transferred  to  a  metallic  collector  probe.   Since  the  col- 
lector is  also  negatively  charged,  work  is  accomplished 
and  some  of  the  kinetic  energy  of  the  gas  has  been  changed 
directly  into  electrical  energy.   These  charges  flow 
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through  an  outside  circuit  and  return  to  the  needle.   If 
the  energy  drained  from  the  gas  is  minimal,  then  the  flow 
has  not  been  disturbed  by  the  conversion  process. 

One  of  the  important  parameters  here  is  particle  size. 
If  the  particle  is  too  small,  i.e.,  very  mobile,  the  drift 
velocity  of  the  particles  becomes  large  and  current  loss 
results  as  exemplified  in  the  current  density  equation  for 
EGD, 

j  =  p^  (U^  -  yE)  (1) 

where  yE  is  the  drift  velocity.   By  rearranging  the  equa- 
tion in  its  scalar,  one -dimensional  form  one  obtains 

j  =  p^  (1  -  HE  3  u^  C2) 

00 

and   slip    can  be    defined   as 

s    =    ^    .  (3) 

oo 

For  a  certain  range  of  radii,  as  the  charged  particle  size 
decreases,  the  mobility  constant  increases  (Figure  2)  and 
hence  slip  increases  [3,9]. 

Mobility  is  a  function  of  mass  and  charge  of  the 
particle.   Since  the  ability  of  a  charged  particle  to  fol- 
low the  turbulent  eddies  in  the  cylinder  wake  is  a  function 
of  its  mass,  an  increase  in  size  results  in  lower  frequency 
response  of  the  charged  particle.   This  can  be  seen  by  con- 
sidering the  effects  of  inertia  in  a  moving  particle.   A 
particle  with  less  mass  can  accelerate,  decelerate,  and 
change  direction  much  faster  than  a  particle  of  greater 
mass  exposed  to  the  same  forces. 
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One  of  the  problems  in  previous  investigation  [6] 
has  been  the  generation  of  steam  droplets  of  correct  size 
for  high-frequency  coupling  and  correct  charge  for  low 

mobility.   Injection  of  particles  ranging  between  10 

-  7 
to  10   meters  charged  in  a  corona  field  results  in  a 

maximum  drift  velocity  of  one  to  ten  meters  per  second 
but  allows  particles  to  fluctuate  with  frequencies  up  to 
50  kHz  [4].   Previous  work  showed  a  high  frequency  drop 
for  the  EGD  probe  in  comparison  with  hot-wire  anemometer 
measurements  [6].   An  analysis  indicated  that  the  particle 
radius  was  on  the  order  of  five  microns  [4] .   Therefore  it 
seems  that  the  size  of  the  steam  droplets  can  be  reduced 
by  decreasing  the  condensation  which  should,  in  turn,  re- 
sult in  an  increase  in  EGD  probe  bandwidth.   Previous  work 
showed  that  "condensation"  was  decreased  by  lowering  the 
steam  pressure  to  below  8.3  psig  [4].   Wetness  could  be 
decreased  by  superheating  the  steam  as  it  travels  from 
the  steam  generator  to  the  cylinder  exit  nozzle.   In  ad- 
dition, Reference  13  (unpublished)  and  some  of  our  pre- 
liminary work  indicated  that  certain  ranges  of  relative 
humidity  of  air  in  the  channel  could  enhance  particle  si^e 
after  injection  into  the  air  flow.   Thus,  high-frequency 
response  could  also  be  enhanced  by  controlling  relative 
humidity. 

One  last  area  concerns  the  standard  of  comparison. 
The  hot-wire  anemometer  measures  the  turbulence  in  a 
single-phase  flow  while  the  EDG  probe  measures  a  two-phase 
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flow.   The  frequency  spectrum  of  these  two  flows  may  dif- 
fer since  the  droplets  most  likely  enhance  the  low  fre- 
quency spectrum. 
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III.   EQUIPMENT 

The  heart  of  the  apparatus  is  the  small  plexiglass 
wind  tunnel  test  section.   The  inside  dimensions  are  two 
by  four  inches  in  cross  section  and  ten  inches  in  length. 
The  test  section  contains  a  steam  injector  unit  which  is 
fed  by  a  heated  steam  line  from  a  steam  generator  [10]. 
The  aerosol  injector  nozzle  consists  of  a  stainless  steel 
ring  embedded  in  the  cylinder  wall  which  allows  the  in- 
jection of  condensing  steam  into  the  turbulent  wake  of 
the  cylinder  (Figure  3) . 

The  steel  ring  also  serves  as  the  corona  ring  posi- 
tive electrode  of  the  aerosol  ionizer.   The  negative 
electrode  consists  of  a  steel  needle  positioned  on  the 
corona  ring  centerline  but  slightly  recessed.   A  Sorensen 
High-Voltage  DC  Power  Supply  provides  the  high  potential 
between  the  ring  and  needle  necessary  to  negatively  charge 
the  steam  particles  as  they  pass  through  the  throat  of 
the  nozzle  (at  sonic  velocity)  into  the  turbulent  cylin- 
der wake.   The  potential  difference  is  monitored  by  a 
Sensitive  Research  electrostatic  voltmeter.   A  few  milli- 
meters downstream  from  the  cylinder  nozzle  a  third  elec- 
trode, consisting  of  a  stainless  steel  probe  with  a  rounded 
tip  (Figure  4) ,  collects  the  current  and  thereby  detects 
the  frequency  spectrum  of  the  turbulent  wake  via  the  infor- 
mation carried  by  the  charged  particles. 

The  air  is  supplied  by  a  Carrier  three-stage  centri- 
fugal compressor  with  a  maximum  flow  rate  of  4000  cubic 
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feet  per  minute  and  a  maximum  pressure  ratio  of  two.   This 
air  is  fed  through  a  cooling  bank  into  the  main  plenum 
which  is  attached  to  the  plexiglass  test  section.   The  valve 
on  the  cooling  bank  which  controls  the  amount  of  air  being 
cooled  also  enables  control  of  the  temperature  and  hence 
the  relative  humidity  of  the  air  (Figure  5) .   The  air  passes 
through  the  test  section  and  around  the  cylinder,  then 
mixes  with  the  steam  jet.   The  particles  are  convected  to 
the  collector  probe,  where  they  are  neutralized,  and  the 
water-air  mixture  flows  out  of  the  test  section  into  the 
atmosphere . 

The  current  from  the  collector  probe  is  fed  into  a 
black  box  which  houses  the  circuitry  for  the  probe  (Figure 
6) .   The  three  International  Rectifier  diodes  prevent  re- 
verse current  flow  [11].   A  Simpson  microammeter  measures 
the  collector  current.   The  capacitor  and  the  neon  bulb 
safety  switch  protect  the  frequency  analyzer.   The  capaci- 
tance is  selected  large  enough  to  ensure  a  phase  shift  of 
less  than  one  degree  in  the  frequency  spectrum  from  25 
hertz  to  20  kilohertz  and  to  accomodate  the  high  voltages 
which  are  present  (Appendix  A) . 

The  pre-amplif ier  (Figure  7)  is  a  self-powered  unit 
consisting  of  a  UA741C  operational  amplifier,  resistors, 
and  a  battery.   It  provides  a  gain  of  ten  over  a  frequency 
range  of  100  kilohertz.   A  shunt  resistance  of  39,000  ohms 
is  placed  across  the  input  in  order  to  avoid  saturating 
the  operational  amplifier  in  the  preamplifier,  i.e.,  one 
third  of  the  signal  is  bled  off  by  this  shunt. 
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The  frequency  analyzer  is  a  General  Radio  Company  Type 
1921  Real  Time  Analyzer  (Figure  8) .   This  consists  o£  a 
1925  Multifilter  and  a  1926  Multichannel  RMS  Detector.   The 
signal  (Figure  9)  from  the  preamplifier  is  put  into  the 
multifilter,  separated  according  to  frequency,  put  into  the 
RMS  detector  where  the  data  are  converted  to  binary  digital 
form  and  the  rms  level  computed  by  true  linear  integration. 
The  computed  data  are  then  displayed  digitally  on  the  rms 
detector  nixie-tubes  and  the  analog  form  of  the  data  is 
projected  graphically  on  a  Hewett -Packard  X-Y  display. 
Either  a  Tektronix  Type  RM  15  or  a  Type  551  Dual  Beam 
oscilloscope  is  used  to  monitor  the  signal  input  to  the 
frequency  analyzer. 

The  hot-wire  data  of  this  experiment  were  collected 
with  a  Thermo  Systems,  Incorporated,  Model  1050  and  1051-D, 
constant  temperature  hot-wire  anemometer  (Figure  10) .   The 
frequency  response  of  the  TSI  equipment  was  superior  to 
the  frequency  response  of  two  other  constant  temperature 
hot-wire  anemometers  which  were  also  tried,  i.e.,  a  Flow 
Corporation  Constant  Temperature  Anemometer,  Model  900, 
and  a  Security  Associates,  Model  100,  hot-wire  anemometer. 
In  addition,  the  TSI  unit  had  a  built-in  signal  generator 
that  enabled  simplified  procedures  for  frequency  calibra- 
tion.  A  Kronhite  Model  3202R  bandpass  filter  was  used  in 
conjunction  with  the  TSI  unit  to  filter  out  the  noise  from 
outside  the  test  frequency  spectrum,  which  was  25  to  20,000 
hertz . 


other  equipment  included  a  psychrometer  for  measuring 
test  section  air  humidity,  a  Simpson  microammeter  for 
measuring  needle  current  and  a  voltage  regulator  for 
regulating  line  voltage. 
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IV.   PROCEDURES 

The  experimental  procedure  developed  from  an  initial 
familiarization  with  the  operation  of  the  apparatus  and  a 
desire  to  start  with  a  verification  of  previous  work.   This 
produced  results  similar  to  those  reported  in  Reference  6. 

The  succeeding  runs  involved  a  change  from  a  0.5 
microfarad,  2000  Vdc  capacitor  to  a  10  microfarad,  2000 
Vdc  capacitor  in  order  to  eliminate  a  phase  shift  in  the 
portion  of  the  spectrum  below  250  hertz  (Appendix  A) .   A 
closer  correspondence  in  the  lower  half  of  the  frequency 
spectrum  now  existed,  between  the  EGD  turbulence  intensity 
(Tp)  and  the  hot-wire  turbulence  intensity  (T,  ) ,  than 
measured  in  previous  work  [6]. 

Asbestos  insulation  was  wrapped  around  the  bare  por- 
tion of  copper  tubing  immediately  below  the  cylinder.   This 
prevented  early  cooling,  and  subsequent  condensation,  from 
forming  particles  of  too  large  a  mass,  which  translates  to 
a  great  inertia  and  a  low  mobility.   During  the  test  runs 
to  check  these  changes,  attention  was  given  primarily  to 
tape  heat  (this  tape  heats  the  line  between  the  steam 
generator  and  the  cylinder  in  the  air  stream),  steam  pres- 
sure, and  applied  voltage.   The  temperature  range  was  be- 
tween 230°F  and  288°F,  the  pressure  between  7.9  to  14.1 
psig,  and  the  voltage  between  215  to  3.1  kilovolts. 

Since  much  of  the  data  taken  on  different  days  showed 
little  correlation  even  when  the  above  parameters  remained 
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unchanged,  it  was  necessary  to  investigate  other  parameters. 
Closer  observation  of  the  relative  humidity,  use  of  a  more 
accurate  device  for  measurement  of  humidity,  and  measurement 
of  humidity  in  the  test  section  during  operation  revealed 
that  the  best  results  occurred  on  days  when  the  relative 
humidity  in  the  channel  was  less  than  461. 

For  operation  at  25^  relative  humidity,  test  section 
temperature  ranged  between  95°F  to  125°F.   Difficulty  was 
encountered  in  getting  values  of  collector  current  and 
total  current  greater  than  zero  when  the  relative  humidity 
was  251.   Even  at  30%,  I^  and  I^  tended  to  be  erratic,  the 
microammeters  registering  zero  approximately  half  of  the 
time.   Good  results  occurred  at  32%  and  this  value  was 
chosen  for  subsequent  testing. 

The  next  portion  of  the  experiment  was  devoted  to  a 
thorough  testing  of  the  equipment  after  it  was  discovered 
that  the  internal  calibration  voltage  was  occasionally  er- 
ratic.  The  General  Radio  equipment  tested  in  accordance 
with  manual  specifications.   Noise  in  the  output  at  the 
low  frequencies  (most  prevalent  around  60  hertz)  was  re- 
duced by  grounding  all  equipment  used  in  the  experiment. 
A  line  voltage  regulator  was  installed  to  decrease  fluc- 
tuations in  line  voltage  which  might  affect  the  equipment 
and  results. 

A  method  of  comparing  existing  hot-wire  results  with 
EGD  measurements  was  devised  using  the  multifilter  unit. 
By  subtracting  out  the  values  (in  decibels)  of  the  turbulence 
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intensity  as  measured  by  a  hot-wire  throughout  the  frequency 
spectrum  on  each  o£  the  multifilter  channels,  any  deviation 
from  a  straight  line  on  the  X-Y  display  gave  an  immediate 
indication  of  the  portion  of  the  spectrum  with  which  the 
EGD  measurements  did  not  favorably  compare  and  the  amount 
by  which  the  intensity  differed  from  the  hot-wire  data. 
Because  it  appeared  that  some  values  of  the  multifilter 
settings  would  subtract  out  a  greater  value  of  turbulence 
intensity  than  current  EGD  measurements  were  producing, 
the  preamplifier  was  incorporated. 

The  hot-wire  anemometer  portion  of  this  experiment 
was  conducted  to  obtain  measurements  of  turbulence  in- 
tensity, T,  ,  at  a  nozzle  pressure  of  eight  psig  and  to 
verify  previous  data  taken  at  13  psig.   In  order  to  simu- 
late the  jet  effects  but  avoid  the  two-phase  flow  in  which 
the  hot-wire  calibration  is  not  useful  and  also  to  avoid 
breaking  the  hot-wires  by  droplet  impact,  the  steam  unit 
was  disconnected  and  air  at  the  same  pressure  was  injected 
into  the  cylinder  and  out  of  the  nozzle  at  8  and  13  psig 
[10].   A  limited  amount  of  testing  (which  occurred  before 
the  failure  of  a  field-effect  transistor)  with  the  Securitv 
Associates  Hot-Wire  Anemometer  produced  results  which  agreed 
with  those  of  reference  [6]. 

A  Thermo  Systems,  Incorporated  (TSI)  constant  temper- 
ature anemometer  provided  a  much  greater  intensity  at  the 
high  frequency  end  of  the  spectrum  than  had  been  previously 
measured  (Figures  11  and  12) ,  especially  when  nozzle  air 
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was  on.   Figure  11  was  the  result  obtained  using  a  DISA 
55F31  probe.   All  succeeding  runs  were  made  with  probes 
used  in  previous  work  which  are  locally  manufactured. 
Subsequent  runs  with  a  bandpass  filter  did  not  improve 
the  results.   It  was  later  discovered,  however,  that  a 
signal  greater  than  0.5  volts  was  possibly  injected  into 
the  General  Radio  1560-P40  preamplifier,  damaging  the 
field-effect  transistor.   Even  though  hot-wire  measurements 
were  made  bypassing  this  preamplifier,  it  did  obtain  power 
from  the  frequency  analyzer  and  was  left  connected  to  the 
analyzer.   It  is  not  known  whether  this  could  have  caused 
a  stray  signal  to  be  injected  into  the  analyzer  causing 
this  unexpected,  high  frequency  "tail."  Another  measure- 
ment with  the  TSI  unit  was  attempted  after  the  above  ex- 
perience and  culminated  in  failure  to  the  TSI  unit  before 
any  results  could  be  obtained. 

The  final  portion  of  this  experiment  produced  the 
EGD  results  reported  herein.   A  smaller  stainless  steel 
probe  (0.062  inch  diameter)  replaced  the  0.125  inch  steel 
probe  used  initially  [6].   It  was  felt  that  a  smaller  probe 
would  sample  a  more  localized  area  and  not  tend  to  average 
the  turbulent  eddies  present  in  a  larger  volume  of  space. 
A  new  preamplifier  (Figure  7)  replaced  the  General  Radio 
1560-40.   A  shunt  resistance  was  incorporated  to  avoid 
over-driving  the  operational  amplifier.   Channel  humidity 
was  held  at  321,  as  measured  at  the  exit  of  the  test  sec- 
tion, eight  centimeters  aft  of  the  cylinder  nozzle.   Tape 
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temperatures  were  varied  between  230°F  and  295°F.   Generator 
collector  current  was  varied  between  zero  and  10.2  microam- 
peres and  applied  voltage  between  zero  and  4.43  kilovolts. 
Collector  currents  less  than  0.2  microamperes  were  not  mea- 
sured. 

The  frequency  analyzer  converted  the  AC  electrical 
input  from  AC  voltage  to  the  appropriate  value  in  decibels. 
The  measure  of  turbulence  intensity,  T,  at  each  frequency 
sampled,  was  expressed  in  dB  and  read  directly  from  the 
X-Y  display  or  the  nixie  tubes  (Appendix  B) .   The  pro- 
cedure developed  for  comparing  the  EGD  data  with  the  hot- 
wire results  is  as  follows:   Tp  represents  the  value  of 
turbulence  intensity  measured  with  the  EGD  probe  as  given 
by  the  frequency  analyzer  and  expressed  in  decibels.   T, 
is  likewise  the  value  of  turbulence  intensity  measured  with 

the  hot-wire  anemometer.   If  the  difference,  AT  =  T,  -  T  t^  , 

'       h    E ' 

is  the  same  (±  1  dB  for  error)  for  each  frequency,  then  the 
difference  results  solely  from  the  difference  in  amplifica- 
tion of  the  two  signals  and  the  two  signals  are  considered 
to  be  identical.   In  this  context  T  =  log, ^  T'(f)  (Appendix 
C). 

The  procedure  developed  for  determining  the  high  fre- 
quency response  involved  a  comparison  of  the  difference  in 
intensity  between  the  most  turbulent  and  the  least  turbulent 

frequencies.   T    is  the  maximum  turbulence  intensity, 
^  max 

usually  found  at  either  the  Strouhal  or  twice  the  Strouhal 

frequency.   T  .   is  the  minimum  turbulence  intensity,  in 
^    ^  mm 
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decibels,  normally  found  'at  the  highest  frequency.   AT„T7n  ~ 

T    -  T  .    is  a  measure  of  the  difference  in  intensity 
max    mm  ' 

of  the  turbulence  between  the  two  frequencies.   If  this 
difference  is  small,  then  a  good  relative  high  frequency 
response  is  indicated.   ^T„pn  will  always  have  a  value 
greater  than  zero  since,  theoretically,  turbulent  vertices 
at  high  frequencies  dissociate  rapidly  to  lower  frequency 
vertices  [7] .   The  typical  intensity  versus  frequency 
curve  for  a  free  jet  exhibits  a  gradual  decrease  in  in- 
tensity as  the  frequency  increases  [7] .   The  best  (smallest) 

value  of  AT,,T-Ti  obtained  for  both  EGD  and  hot-wire  data  was 
HFR 

12  dB. 
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V.   RESULTS  AND  DISCUSSION 

Good  results  were  consistently  obtained  when  little 
or  no  tape  heat  was  applied  and  when  the  collector  current 
was  very  carefully  controlled  (Table  I) .   Figure  13  shows 
a  corresponding  drop  in  high  frequency  response  as  temper- 
ature is  increased.   In  order  to  examine  and  compare  data, 
a  plot  of  pressure  versus  density  for  various  values  of 
temperature  (Figure  14)  was  constructed  from  steam  tables 
[1,8].   This  graph  was  used  to  convert  values  of  pressure 
and  temperature  into  corresponding  values  of  density.   A 
check  of  the  steam  pressure  and  temperature  values  imposed 
during  the  experiment,  on  this  plot,  reveals  that  best  re- 
sults were  obtained  when  points  lay  in  the  saturated  or 
supersaturated  regions. 

Two  methods  of  comparing  EGD  and  hot-wire  data  were 
used,  both  using  the  data  shown  in  Figure  11  as  the  stan- 
dard of  comparison.   Figures  15  through  17  show  typical 
results  at  12  to  13  psig  with  varying  temperatures  and 
voltages.   With  the  tape  heat  turned  off,  the  temperature 
of  the  steam  at  this  pressure  is  240°F.   Saturation  tem- 
peratures of  the  steam  for  these  pressures  are  244°F  and 
246°F  respectively.   Figures  18  through  20  show  typical 
results  at  8  psig  at  various  temperatures  and  voltages. 
At  this  pressure  the  saturation  temperature  is  235''F  and 
the  steam  temperature  with  the  tape  heat  off  is  230°F. 
These  results  are  compared  with  the  best  hot-wire  anemometer 
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data  available.   These  data  were  chosen  as  the  standard  of 
comparison  since  these  were  the  first  data  collected  with 
the  TSI  unit  and  are  expected  to  be  free  from  any  effects 
due  to  subsequent  malfunctioning  electronics.   Since  all 
data  are  already  converted  to  decibels,  to  ratio  (see  Ap- 
pendix C)  the  desired  points  one  simply  takes  the  dif- 
ference of  the  logarithm  of  the  plotted  points  at  each 
frequency,  measured  from  some  arbitrary  reference  line,  to 
the  points  to  be  matched.   Twice  the  Strouhal  frequency 
was  chosen  as  a  reference  point  since  it  is  the  most  dom- 
inant point  on  the  EGD  and  hot-wire  data  collected  and 
also  since  the  frequencies  just  beyond  this  point  are  most 
responsive  to  changes  in  parameters. 

The  second  method  of  comparing  EGD  data  with  hot- 
wire data  was  by  inserting  the  hot-wire  data  into  the  multi- 
filter  as  explained  earlier.   Figures  21  and  22  show  a 
comparison  of  typical  results  displayed  during  testing. 

Figure  23  shows  the  effect  of  the  density  of  the 
steam  on  the  current  efficiency  (1^/1^^)  for  the  two-phase 
fluid.   In  addition  to  the  sensitivity  of  the  current 
efficiency  to  a  small  change  in  voltage,  the  effect  of  an 
increase  in  density  is  to  require  a  higher  voltage  in 
order  to  obtain  the  same  current  efficiency.   This  shift- 
ing of  the  curves  to  the  right  is  similar  to  that  reported 
by  Whitby  [12]  for  the  effect  of  pressure  (temperature  was 
held  constant)  in  a  single-phase  fluid  (air) . 
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Comparisons  of  Table  I  with  Figure  13  and  Figures  15 
through  22  reveal  that  the  best  high  frequency  response, 
i.e.,  when  AT^^pj^  is  least,  generally  occurs  when  the  steam 
temperature  is  at  or  below  saturation  temperature  for  a 
given  pressure.   Figure  24  is  a  plot  of  pressure  versus 
high  frequency  response  for  three  ranges  of  temperature 
to  saturation  temperature  ratios.   The  relative  humidity 
and  the  collector  probe  distance  remained  constant.   The 
corona  voltage  was  held  constant  at  three  kilovolts  with 
exceptions  shown  beside  the  data  points.   The  best  high 
frequency  response  (lowest  AT^p^)  occurs  at  the  super- 
saturation  and  saturation  temperatures  with  the  best  re- 
sults in  the  eight  psig  pressure  range. 

Figure  25  shows  that  a  definite  optimum  voltage  exists 
for  operation  at  good  high-frequency  response.   All  other 
variables  are  held  constant  but  a  slight  decrease  in  the' 
relative  humidity  from  341  to  32%  and  a  slight  increase 
in  pressure  from  8.1  to  8.2  psig  during  the  experiment 
may  have  had  an  effect  on  the  optimum  collector  current 
and  on  the  current  efficiency.   The  accompanying  changes 
in  collector  current  and  current  efficiency  are  shown. 
After  a  "triggering"  voltage  is  attained,  collector  cur- 
rent increases  linearly  until  finally  a  plateau  is  reached. 
The  graph  shows  that  the  optimum  current  efficiency  does 
not  occur  at  the  optimum  voltage,  although  some  previous 
observations  indicated  that  good  results  occurred  at  a 
current  efficiency  of  0.2  to  0.5  and  at  a  higher  collector 
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current.   Figure  26  shows  the  effect  of  current  efficiency 
for  three  different  temperatures  with  other  parameters 
held  constant. 
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VI.   CONCLUSIONS 

It  is  concluded  that  an  electrogasdynamic  probe  can 
indeed  by  successfully  used  to  measure  intense  turbulent 
fields.   Results  comparable  to  or  better  than  some  of  the 
existing  hot-wire  anemonmeters  have  been  obtained. 

The  following  specific  conclusions  are  drawn: 

1.  A  controlled  relative  humidity  of  32%  resulted  in 
a  significant  improvement  in  the  measurement  of  turbulence 
intensity  in  the  high  frequency  range  of  the  spectrum. 

2.  Little  or  no  tape  heat  gave  optimum  results  (Fig- 
ure 22).   This  parameter  is  a  highly  critical  one.   A 
change  of  1°F  can,  at  times,  make  a  difference  of  measuring 
a  highly  intense  turbulence  or  not  detecting  any  turbulence 
at  20  kHz.   Sometimes  saturation  temperature  provides  op- 
timum results.   At  other  times  no  tape  heat  (5°F  below 
saturation  temperature)  provides  optimum  results.   Even  in 
this  last  case,  an  increase  in  tape  heat  of  one  degree  re- 
sults in  a  marked  steadiness  of  total  current  and  collector 
current.   An  increase  of  an  additional  1°F  would  steady 
the  current  even  more  but,  for  some  runs,  would  also  res:  It 
in  a  significant  high  frequency  detection  loss.   This  fact 
is  contrary  to  expectation,  (earlier  investigation  indicated 
that  the  steam  droplets  seemed  too  large  already) . 

3.  A  ratio  of  collector  current  to  total  current 
greater  than  0.20  seemed  to  provide  best  results  (Figure  26). 

4.  A  marked  improvement  could  sometimes  be  seen  dur- 
ing transient  conditions,  especially  after  a  change  in  voltage 
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The  time  required  to  reach  steady  state  conditions  is 
generally  30  seconds  to  one  minute.   However,  in  a  few 
instances,  transients  were  noted  even  after  five  minutes, 
indicating  that  steady  state  conditions  are  important. 

5.  An  optimum  applied  voltage  range  exists,  below 
which  no  collector  current  is  produced  and  above  which 

no  high  frequency  turbulence  can  be  measured.   This  latter 
phenonemon  may  be  due  to  an  increasing  charge  built  up 
per  particle,  resulting  in  an  increased  drift  velocity 
back  upstream  to  the  attractor  ring  (Figure  25) . 

6.  The  effect  of  an  increase  in  steam  pressure  is 
to  shift  the  current  efficiency  curves  to  the  right 
(Figure  23) .   This  is  comparable  to  conclusions  drawn  by 
Whitby  [10]  for  single-phase  flow.   Saturation  effects  of 
the  steam  were  lost  in  converting  pressure  and  temperature 
to  density,  however,  and  are  not  accounted  for  in  this  fig- 
ure.  These  also  seem  to  play  an  important  role. 

7.  Eight  psig  seems  better  than  thirteen  psig.   This 
is  in  agreem.ent  with  our  previous  arguments  for  smaller 
mass  particles. 

8.  Total  collected  current  seemed  to  have  an  effect. 
For  best  results  a  collector  current  of  one-half  microam- 
pere or  greater  was  required  (Table  I) . 
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VII.   RECOMMENDATIONS 

Several  areas  need  to  be  further  explored  in  order 
to  refine  the  electrogasdynamic  probe  and  enhance  future 
applications  of  this  device.   The  areas  are  the  follow- 
ing: 

1.  Measurement  of  the  turbulence  intensity  with 
state-of-the-art  hot-wire  anemometer  equipment  to  deter- 
mine an  accurate  and  consistent  turbulence  intensity  fre- 
quency spectrum. 

2.  An  investigation  into  the  effect  of  higher  ap- 
plied voltages  on  the  production  of  charged  particles  and 
a  better  high  voltage  supply,  i.e.,  regulated  to  avoid 

±  60  volt  fluctuations  at  3.0  kilovolts  (RMS  ripple  factor 
is  2.01)  . 

3.  Better  instrumentation  for  the  control  of  steam 
pressure  and  temperature,  i.e.,  more  accurate  monitoring 
of  these  data.   Temperature  control  should  be  within  at 
least  ±  0.5°F,  preferably  ±  0.2°F  and  pressure  control  to 
at  least  ±  0.05  psi.   This  would  allow  accurate  deter- 
mination of  steady  state  conditions  and  insure  that  the 
critical  parameters  are  indeed  constant.   Also,  continuous 
monitoring  of  air  humidity  upstream  of  the  cylinder  is 
desirable  to  eliminate  interruption  of  the  experiment  for 
humidity  measurement  and  to  preclude  any  unknown  changes 
from  occurring  between  measurements. 
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4.  A  study  to  determine  whether  there  is  any  signi- 
ficant difference  between  the  turbulent  field  of  a  single- 
phase  fluid  and  that  of  the  fluid  when  a  second  phase  is 
introduced.   It  is  possible  that  the  two-phase  flow  may 
shift  the  spectrum  towards  lower  frequencies  by  virtue 

of  the  large,  massive  particles  which  are  present. 

5.  Elimination  of  the  temperature  difference  be- 
tween the  free  stream  and  the  nozzle  air.   In  addition  to 
the  lack  of  two-phase  flow,  using  the  hot-wire  as  a  ref- 
erence introduces  another  discrepancy  -  the  presence  of 
two  temperatures  in  the  fluid.   Since  the  hot-wire  mea- 
sures temperature,  any  sensing  of  intermittent  tempera- 
tures due  to  two  different  temperatures  rather  than  to  a 
difference  in  turbulence  intensity  could  lead  to  a  signi- 
ficant distortion  of  the  frequency  spectrum  curve. 

6.  A  theoretical  study  of  the  nature  and  frequency 
range  of  the  turbulence  in  the  wake  of  the  cylinder  and 
nozzle  to  determine  if,  indeed,  the  hot-wire  anemometer 
is  measuring  all  effects  present. 

7.  Further  investigation  of  the  effect  of  small 
changes  in  temperature  and  pressure  on  the  formation  of 
droplets  and  the  production  of  collector  current  at  various 
frequency  ranges  in  the  turbulence  spectrum.   An  accurate 
determination  of  the  size  of  the  droplets  as  a  function 

of  the  various  parameters,  especially  temperature  and 
ppressure,  at  a  given  humidity,  will  greatly  aid  in  deter- 
minimg  the  proper  environment  for  the  future  use  of  the 
EGD  probe. 
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FIGURE   4. 
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FIGURE     6. 
ELECTRICAL    SCHEMATIC     OF    COLLECTOR 
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FIGURE    7. 
SCHEMATIC    OF    PREAMPLIFIER 


o  OUTPUT 
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X-Y   DISPLAY,   RMS   DETECTOR    AND    MULTIFILTER 
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MULTIFILTER    UNIT 


FIGURE    8 
FREQUENCY     ANALYZER 
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TSI    HOT-WIRE  ANEMOMETER    AND    FILTER 


Apparatus  Ready    for  Testing 

FIGURE    10 
CONFIGURATION     FOR    HOT-WIRE   TESTING 
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MAIN    AIR     ONLY   (IBo  rps) 


MAIN    AIR(l85fps)    NOZZLE   AIR  (8psig) 


FIGURE    II. 


X-Y  DISPLAY    FROM    HOT-WIRE   ANEMOMETER 
L=  10  mm   (1) 


43 


MAIN    AIR   ONLY     (185   fps) 


MAIN    AIR    (185  fps)     NOZZLE    AIR  (8psig) 

FIGURE   12.    X-Y  DISPLAY   FROM    HOT-WIRE    ANEMOMETER 

L  =  5  mm    IQ_) 
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=  0.9   uA,    MAIN    AIR(l85fps),   STEAM   (8.0  psig,  245  °F) 


Reference    — 


lr  =  0.9aA,    MAIN    AIR  (185  fps) ,  STEAM   (8.0  psig,  231  °F) 


FIGURE    13.     X-Y    DISPLAY    FROM     EGD    PROBE 

L  =  5   mm  (  C^.) 
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FIGURE     2! 
COMPARISON    OF    EGD    DISPLAYS 
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COMPARISON    OF  EGD     DISPLAYS 


54 


0,8 


.7 


m 


Q 


8   psig  — «-' 


e 


ffl" 


el 


e 


1 
I 


s 


o 


0 


SMALL    PROBE    DATA 

RH  =  32  % 
L    =    5     mm 


I 
\ 

ttM? 


e 

A 
O 

D 
O 


0.055 
0.056 

0.064 
0.065 
0.066 
0.067 


Ib^/ft^ 


8    psig 


12-13  psig 


^ 


^ 


*• —  12"  13   psig 


a 


\0 


o 
o 


A 


0.068 


A 


aA 


0 


1.0 


4.0 


± 


5.0 


2.0  3.0 

Vc    (kV) 

FIGURE    23.      EFFECT     OF    DENSITY     ON     CURRENT    EFFICIENCY 


55 


o 

V) 


<D  O  O 

6  -  -■ 

I      I     I 

^  cn 
d  O- 


DOO 


a>' 


o  o 
o. 
c 
$  o 

-I 

«/)   <l> 

c  <w 


0^ 

(VJ 

to 
II 

I 


E> 

o 

II    II 


> 


(vi 

o 


< 


O    rs 

CD 
□ 


D 


D     D 


n 


D 


<3 


00 


fO 

<] 


<]03  O 


03     □ 

O 


o 

ro 
o 

LiJ 

(r. 

cvi 

CO 
CO 
LU 

a. 

O 

Q 

2 

< 
UJ 

a: 

O 

d 

1- 
< 

o> 

q: 

Ui 

LJ 

n 

n 

hi 

o 

Q. 

1- 

0) 

o 

GO 


o 
^-■ 


i 


o 


o 

ro 


O 


CL 
O 


CVJ 

UJ 

ID 

Ll 


XJ 


56 


_C>j_H 


d 

6 

ro 
d 

CVJ 

d 

d 

o 

1 

1 

1 

1 

1 

•  1 

o 

a> 

Q.0     H-    C  O^ 


o 

ro 


_?      i 


o 
cvi 


o 


rsj 

00 


fO 
,   Oin       I 

oocM  —      ro 
II  (I  II   II    It 

Q.H3-JQ: 


IT) 


00 


CJ 


'^ 

a: 

OJ 

u. 

^-^ 

I 

■ 

m 

1— 

•u 

*" 

< 

CJ 


O 

I" 


lO 


o 


o 

ro 


> 


>" 


LU 
O 


LlI 
(3 
< 

H 
_I 

> 


H 
Ol 
O 


lO 
CVJ 

LU 
01 

3 
C£ 

li. 


•-      oQ-i — 235       I— OX      ii.Q:iiJ03LiJZo>       trujc/jQ-Ozcrtuj 


57 


u. 

o 

O  — lO 

ro  fO  ro 

CO  CVJOJ 


I 
I 
I 


0> 
(A 


0> 

JC 


o 


o  to 
0)  X3 


/  tv- 


o> 


o 

0) 

CO 
(M  ,.  »C^  O^  O 
M^  lO  CO     •  — c 
fO        _  CVJ  CO 

II    II    II    II    II 

X    ,   8 
en- 


>Q- 


J. 


O 


O 
ro 


lO 

O 

CVJ 

a: 

CJ 

u. 

^— » 

.1 

CD 

II- 

XJ 

< 

«o 


f^ 

o 

6 

2 

Ix. 
< 

CD 

d 

to 

< 

1x1 
(/) 

2 

lO 

O 

6 

Y     RESP 
NCY 

^ 

, 

O  UJ 

6 

CD 
CM 

go 

j- 

3  Li- 

v^ 

liJ 

O  Lx 

ti) 

q: 

LU  Ld 

3 

cr 

ro 

CD 

^h- 

6 

u. 

N     HIGH 
CURREN 

CVJ 

o 

6 

LU 
(T 

3 

< 

*■"■ 

UJ 

d 
o 

Q_ 

UJ 

U_ 
O 

CJ 
UJ 

UJ 


58 


APPENDIX  A:   PHASE  SHIFT  CALCULATIONS 

The  phase  shift  is  calculated  below  in  the  order  in 
which  equipment  changes  were  made.   The  relation  governing 
phase  shift  is : 

R 


1. 


tan    (|)   =   -^ — 
1 


c)) 


ooCR 
1 


27TfCR 
C   =    0.5  microfarads,    R  =    100,000   ohms. 
At    f  =    25    Hertz : 

-1  1 


(j)   =    tan 


2    (25)  (0.5    X   10"^)  (10^) 


=  tan  ■"  (0.127) 
=  7.3° 
At  f  =  20,000  Hertz: 

({)  =  tan"-"-  (0.00016) 
=  0° 
C  =  10  microfarads,  R  =  100,000  ohms. 
At  f  =  25  Hertz : 

(})  =  tan"-'-  (0.0064) 
=  0° 
C  =  10  microfarads,  R  =  825,000  ohms. 
At  f  =  25  Hertz: 

4)  =  0° 


_1 

(joC 
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APPENDIX  B:   DETERMINATION  OF  V 


RMS 


The  signals  measured  in  this  investigation  are  ex- 
pressed in  decibels  for  ease  in  manipulation.   The  decibel 
is  simply  a  ratio  of  some  power  level,  W,  to  a  reference 


power  level,  W  ,  and  defined  as 
o 


Vrms  =  ^0  i°gio  r- 

o 


=  10  lo) 


o 


=  10  log^Q  (V/V^)^ 

=  20  log^Q  (V/V^)  dB. 

A  zero  dB  gain  implies  that  V  =  V  .   If  V  =  lOV  then 
^      ^  o  o 

Vn,4c  =  20  dB.   Likewise,  if  V  =  100,000V   then  V^..^,  =  100 
RMS  '  '     o        RMS 

decibels . 

The  signal  into  the  RMS  detector  is  integrated  over 
a  four-second  period  during  which  each  of  thirty  discrete 
frequencies  are  sampled  simultaneously  1024  times  in  one- 
third  octave  bands.   The  detector  measures  values  of  vol- 
tage with  a  one-millivolt  reference  level.   This  is 
converted  to  decibel  form  using  the  equation 


N  r 


V 


RMS 


=  10  l^SlO   ^  Z 


i=l  L 


V(t.) 


V 


dB 


where 


V(t.)  =  measured  value 
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V  =  reference  level 
o 

N  =  number  of  samples  (binary  number) . 
The  computed  data  are  then  displayed  digitally  on  the  nixie 
tubes  and  are  projected  on  the  X-Y  display  in  analog  form. 

A  scale-factor  adder  adjusts  the  output-level  read- 
ings.  Normally  a  one-volt  RMS  input  produces  an  output- 
level  reading  of  60  dB .   This  reading  is  adjustable  in  one 
dB  steps  to  a  maximum  of  159  dB .   This  adjustment  does  not 
affect  the  analog  output  levels. 

Thus,  in  order  to  determine  V„j,,„ ,  the  nixie  tubes  or 
the  X-Y  display  are  read  as  follows: 
1.    Direct  Reading  from  X-Y  Display. 

r 


T  = 


1"  grid 


-Calibration:  2"  =   \N 


Reference  level 


Reading  from  the  Nixie  Tubes . 


X-Y  display 

Nixie  tubes 

Remarks 

40  dB 

40  dB 

This  reading  is  present 
when  no  signal  is  injected. 

40-80  dB 

40-80  dB 

Actual  values  on  nixies 
are  correct 

80-120  dB 

40-80  dB 

Add  40  dB  to  value  shown 
on  Nixies  to  obtain  actual 
value 

61 


The  frequencies  at  which  measurements  are  made  are 
displayed  both  along  the  x-axis  of  the  X-Y  display  and  on 
the  nixie  tubes  in  the  following  sequence: 


Band 

Frequency 

Band 

Frequency 

Band 

Frequency 

14 

25    Hz 

24 

250   Hz 

34 

2500  Hz 

15 

31.5 

25 

315 

35 

3150 

16 

40 

26 

400 

36 

4000 

17 

50 

27 

500 

37 

5000 

18 

63 

28 

630 

38 

6300 

19 

80 

29 

800 

39 

8000 

20 

100 

30 

1000 

40 

10000 

21 

125 

31 

1250 

41 

12500 

22 

160 

32 

1600 

42 

16000 

23 

200 

33 

2000 

43 

20000 
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APPENDIX  C:   TERMINOLOGY  RELATIONSHIPS 

From  Reference  4,  the  ratio  of  the  EGD  probe  output 
divided  by  the  corresponding  U^m^c ,  and  then  normalized, 
results  in 

I '/An  q 
s        ^~RMS 

where 

I '  =  spectral  probe  current 
A  =  area  of  the  collector 
n  =  number  density  of  charged  particles 
q  =  charge  per  particle 
'^ooRMS  ~  ^^  value  of  the  fluctuating  component  of  carrier 
gas  velocity. 

Let  e  =  reference  voltage  of  one  millivolt.  Dividing  both 
the  numerator  and  the  denominator  by  e  and  taking  the  log- 
arithm of  both  sides  results  in 


where  now 


Tp  =  intensity  of  the  turbulence  as  measured  by  the 


EGD  probe  and  expressed  in  decibels 

intensity  of  the  turbulence  as  measi 

hot-wire  anemometer  and  expressed  in  decibels 


T,  =  intensity  of  the  turbulence  as  measured  by  the 
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